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Abstract 

Wave packet motions of a single electron in harmonic potentials or a magnetic field 
are obtained analytically. The phase of the wave function which depends on both time 
and space is also presented explicitly. The probability density of the electron changes 
its width and central position periodically. These results are visualized using computer 
animation techniques. 


1 Introduction 

We investigate a time evolution of the electron wave packet through analytical methods. The 
time evolutions of restricted initial wave packets were obtained [l]-[3]. Here, we consider a 
general initial wave packet and obtain a classical harmonic oscillation of the center of mass of 
the probability density and an oscillation of its variance. We have also obtained the analytic 
form of the phase of the wave packet. 


2 One-dimensional harmonic potential 

We consider the Schrodinger equation for the one-dimensional harmonic potential 


^ dt^ 2 mdx 2 


hd 2 k 2 
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The stationary solution is 

<p n {x,t) = u n (z)exp(-Mn+ ^)0> 


( 1 ) 

( 2 ) 
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where u> = \fi and «»(*) is expressed using the Hermite polynominal H n (x) 


u n (x) = N n H n (ax) exp(— ^a 2 i 2 ) 


( 3 ) 


with a = and the normalization factor N n (x) = 

Next, we shall expand an initial wave function by these functions and trace its time evolution. 
Hereafter, the unit length oc = 1 is used. Without loss of generality, we choose the initial wave 
packet as 

^ ^ (t z\^ 7 2 

exp( - L i^ L)ex p(-^)’ < 4 > 


ip{x,0) = 


yj 2 tt i 


7T<7 


where 2 is a complex number 2 = z 0 + We shall expand this wave packet in terms of the 
stationary solutions 


#*>0)= £<?»«»(*)• 


( 5 ) 


n=0 


We calculate a expansion coefficient C n , with the help of the generating function of the Hermite 
polynominal and obtain the following expression. 


C n = 


4cm! 


exp(- : 


2 [jl ( 1— 2cr 2 \m/ 2 z \n-2m 


\/27r2 n (l + 2 a 2 ) rv 2(1 + 2a 2 ) 4a 2/ T J^ 0 m!(n-2m)! 

Thus we obtain the time evolution of the wave packet by the following infinite series 


( 6 ) 


rp(x,t) = 


4 a 


\ v'Sfl + 2<r») exp( 
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n = 0 


W~w )“5^ )exp( "T )exp( -T ) 

exp(-ia,()„ (-{^r(l*r)"- 2 ™ 


-r E 


m=0 


l + 2<r 2 / Vl+2<7 2 - 

m!(n — 2m)! 


( 7 ) 


When a 2 equals this summation is evaluated easily. 

Jl z 2 x 2 2 2 

tf'OM) = W ~ ex P(“y ~ ~2 ~ y ) ex P( — 2") ex P(~T ex P(- 2 * w 0 + 2z exp(— iud)). (8) 


The center of mass of the probability density of the wave packet oscillates sinusoidally. On the 
other hand, the variance of the probability density is constant during the motion. 

When the variance a 2 is not we shall eliminate the time dependent phase factor in eq.(7) 
by the following transformations. 


1 - 2cr 2 
1 + 2a 2 


exp(— 2iut) 


1 - 28 
1 + 2.5’ 
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1 - 4<5 2 ' 
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( 10 ) 
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From eqs.(9) and (10) we obtain 


, f 4a 2 + i(l — 4a 4 sin(2w<)) 

6 = So + iti = 2 ((i -f 4cr 4 ) + (1 — 4<j 4 ) cos(2a'<)) ’ 

2 (z 0 cos (wt) + 2zja 2 sin(u><)) + 2t(zj cos(w<) - 2 z 0 a 2 sin(urf)) 
w = w 0 + tw 1 - ( 1 + 4(74) + ( 1 -4a4) COS (2^) ' 

Inserting these values into the expression (7), we have 

~~ ^2 


rp{x } t) = 
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y/2ir{ \ + 2a 2 ) 


exp(- 


y 2 «*2 lUjt 

1 )exp(-y)exp(-— ) 


2(1 + 2a 2 ) 4a 2 
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xT H (*) — V 1+26 1 +-?-*- 

V nlJ 2" L fl m!(n-2m! 

n=0 m=U v / 


( 11 ) 

(12) 


(13) 


Comparing this and the expression (7) at t = 0, we see that this is also the expansion formula 
of a Gaussian wave packet. After a straight forward but lengthy calculation we obtain 


ip(x,t) = exp (17) 


yfio , (*-«>)* / w i \ 

'JyS|«| CXP( « )eXp( 4^o ' 


(14) 


where exp(i7) is a phase factor which depends only on time 

-(z 2 - z 2 ) sin(2wf) + 4z 0 z 1 a 2 (cos(2u;t) - 1) 1 tan (iaQ x 

7 “ 2(cos(2u/<)(l - 4a 4 ) + 1 + 4a 4 ) 2 ^ aiH 2a 2 


(15) 


1 .11 . • 1 n i 1 \J * 2 o ai + *F j V 42 o <j4 +*F 

The center of mass of the wave packet oscillates sinusoidally between — — 5^ — ancl 2^ ■ 

The variance of the probability density changes periodically in the range between and a 2 

( <j 2 > |) or between a 2 and j^y ( a 2 < |). The period of its change is half of that of the 

oscillatory motion of the center of mass[5]. The motion of the probability density function is 

presented in FIG. 1. 



la! <*>l 

FIG. 1. The motion of the probability density function. Here, we choose the 
variance of the probability density of the initial wave packet as (a) bird s-eye 
view, (b) contour line. 
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3 Two-dimensional Harmonic Potential 

Next, we consider the two dimensional Schrodinger equation for an isotropic harmonic potential 

. t d , h 2 V 2 , k 2 (x 2 + y 2 ) 

,h di’ f ' = -^r ,p + ~~2 — *■ (16) 

We choose the initial wave packet as 

ip(x,y, 0) = N exp(—£(x — xo) 2 + ik x o(x~ x o) — r l{y~yo) 2 + ik y0 (y — yo) + \(x — x 0 )(y — yo)), (17) 
where £, 77 and A are complex constants, 

£ = 6 + *’ 6 » V = Vo + , A=A 0 + iA 1) (18) 

which satisfy the following inequalities 

£0 > 0, t?o>0, 4£ 0 t7o — Aq > 0, (19) 


and N is a normalization constant 


/ 4£o*7o — Aq 


Using the same techniques and procedure in the one-dimensional case, we obtain the time 
evolution of the wave packet in terms of an infinite series. 

OO OO 

4>( x,y,t) = £[u»(y) £ Cm,nu m (x)exp(-iu;(m + n+ 1)0], (21) 

n=0 m = 0 

The expansion coefficients C m n are also calculated explicitly. 

For an uncorrelated initial condition 


xp(x, y, 0) = \J v/ ^° 7?0 exp(— £(j - x 0 ) 2 + ik x0 (x - i 0 ) ~ v(y ~ Vo) 2 + ik yQ (y - y 0 )), (22) 


we can evaluate the infinite series 




exp( — ^ ^ 2 _ ' l". )..exp (ik xt x + ik yt y) 


~ tAUl - ~ 

2na xt a yt 4 o* xl 


,.( 1 - 4 <xJ) sin(2w< + 2 j x ){x - x t ) 2 , .(1 - 4 < 7 *)sin( 2 u/f + 27 y ){y - Vt ) 2 ^ 

* 6XP(1 16^5 + 2 16^5 } 

/ k x tXt .k x 0 X 0 .1 ^tan(wt) ^ ^ 

x exp(— t— *— 1 - arcta n ( — ) + id x ) 

, k yt y t .k y oyo .1 . 1 -a \ /oo’i 

x exp(-7— 1— — t - arctanf — ■ - ) + tv v ), (2JJ 
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where 


2 __ (l-4(g + ff))fr 

° x (—16(^o + 6) 2 + 4(£g - £?)) cos 2 ( 7i ) + (1 - 4(£g - ^ 2 ))sin 2 ( 7l ) ’ 

2 (1 -4(T 7 2 + r ? 2 ))77o 

ff * ( — 1 6(?7o + rj{) 2 + 4(j7q - r? 2 )) cos 2 ( 7y ) + (1 - 4(r^ - rfi)) sin 2 ( 7y ) ’ 

\ L , 4& N _ 1 . , 4th , 

7,--wtan( 1 _ 4 ( g + g)) , T»- 2 Mcla "( 1 _ 4W + m 2 ) ). 

2 sin 2 (u4 + 7 i) + 4cr* cos 2 (ujt + 7x ) 2 sin 2 (a;< + 7y ) + 4a* cos 2 (ad 4- 7y ) 

'** = 4^ ’ ^ ” 4cr 2 

Here 9 X and 6 y are time independent phase factors 


e x = 

9 


-( 




(1 - 4cr*)sin(2 7l ) 
4(sin 2 ( 7l ) + 4 <t^ cos 2 ( 7i )) 
(1 ~ 4(7^) sin(2 7y ) 


4(sin 2 ( 7y ) + 4 o\ cos 2 ( 7j/ )) 


+ 6)*o + \ arctan( 

+ 9\)vl + \ arctan( ta ^^ ). 


We obtain the explicit time dependence of the following parameters 


(24) 

(25) 

(26) 

(27) 

(28) 
(29) 


x t = x 0 cos(c A) + k x0 sin(ud) , y t — y 0 cos(a >t) + k y0 sin(od), (30) 

k xt — k x0 cos(uit ) — x 0 sin(ud) , k yt = k yQ cos(u>t) — y 0 sin(od) . (31) 

The trajectory of the center of mass of the probability density function is an elliptic motion 
around the origin with an angular frequency u>. 


4 Uniform magnetic field 

The Schrodinger equation for a single electron in a uniform magnetic field perpendicular to the 
two dimensional flat plane is 

)V, (32) 

ot 2m c 

where the vector potential A in Landau gauge is 

A = (-By, 0) (33) 

We separate a special solution of the wave equation as 

xp(x,y,t) = exf>(ikx)f(k,y,t) (34) 

The wave equation for f(k, y , t) becomes 

,h Jt l{k ' y ' t)= L ( ~h + " a )2)/( *' y ’ (M) 
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where „ 

eB 

a ~ ch 

This is the one dimensional Schrodinger equation for the harmonic potential centered at y = 
kfa. Thus above techniques and procedures can be applied in order to obtain the time evolution 
of the wave packet [6]. We choose the initial wave packet eq.(17). The comlete descriptions 
are presented in the literatutre [6]. The major difference between two dimensional isotropic 
harmonic potential and magnetic field is the period of the change of the variance. The former 
is the half of the latter. This fact is also interpreted by the pass integration technique [6], [8]. 
For the following initial condition 

? = >! , A= j. (37) 

the shape of the contour lines of the probability density function remains circular during the 
motion. 

For the following initial condition 

£ = *7 = ^ , A = ^> ( 38 ) 

the shape of the probability density remains unchanged. 

5 Conclusion 

Using a frame buffer NVS2000 and video recorder BVW-75, we have made CG animations 
which can give us an intuitive understanding of the wave packet motions. 

The potentials are simple but due to the quantum mecanical property the analytic form of 
the wave packet motions are very complicated. 
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